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The process of particle acceleration by left-hand, circularly polarised inertial Alfven waves (IAW) 
in a transversely inhomogeneous plasma is studied using 3D particle-in-cell simulation. A cylindri- 
cal tube with, transverse to the background magnetic field, inhomogeneity scale of the order of ion 
inertial length is considered on which IAWs with frequency 0.3u C i are launched that are allowed to 
develop three wavelength. As a result time-varying parallel electric fields are generated in the den- 
sity gradient regions which accelerate electrons in the parallel to magnetic field direction. Driven 
perpendicular electric field of IAWs also heats ions in the transverse direction. Such numerical 
setup is relevant for solar flaring loops and earth auroral zone. This first, 3D, fully-kinetic simu- 
lation demonstrates electron acceleration efficiency in the density inhomogeneity regions, along the 
magnetic field, of the order of 45% and ion heating, in the transverse to the magnetic field direction, 
of 75%. The latter is a factor of two times higher than the previous 2.5D analogous study and is in 
accordance with solar flare particle acceleration observations. We find that the generated parallel 
electric field is localised in the density inhomogeneity region and rotates in the same direction and 
with the same angular frequency as the initially launched IAW. Our numerical simulations seem 
also to suggest that the "knee" often found in the solar flare electron spectra can alternatively be 
interpreted as the Landau damping (Cerenkov resonance effect) of IAWs due to the wave-particle 
interactions. 

PACS numbers: 52.35.Hr; 52.35.Qz; 52.59.Bi; 52.59.Fn; 52.59.Dk; 41.75.Fr; 52.65.Rr 



I. INTRODUCTION 



Super-thermal particles play an important role in 
many space plasma situations. The relevant two exam- 
ples are: (i) Earths Auroral zone (AZ) that is known to 
host strong field- aligned currents, parallel electric field 
and accelerated particles. Observations essentially show 
two modes of particle acceleration present in AZ: a) Pre- 
cipitating auroral electrons narrowly peaked at specific 
energy, suggestive of a static potential drop in the AZ 
(e.g. Ref.yiJ); b) More recent observations by FAST 
spacecraft (e.g. Ref.[2j]) indicate existence of electrons 
with broad energy and narrow in pitch angle distribution 
that is consistent with the inertial Alfven wave (IAW) 
acceleration, (ii) In solar corona, a significant fraction 
of the energy released during solar flares is converted 
into the energy of accelerated particles jH. The par- 
allel electric field that can accelerate electrons is pro- 
duced when low frequency (u) < u>d, where uj C i — eB/rrii 
is the ion cyclotron frequency) dispersive Alfven wave 
(DAW) has a wavelength perpendicular to the back- 
ground magnetic field comparable to any of the kinetic 
spatial scales such as: ion gyroradius at electron tem- 
perature, p s = \JkBT e /mi/ui C i, ion thermal gyroradius, 
Pi = \J kBTj/mi /uj c i , [J] or to electron inertial length 
A e = c/ujp e [5[. Dispersive Alfven waves are divided into 
Inertial Alfven Waves or Kinetic Alfven Waves (KAW) 
depending on the relation between the plasma /3 and elec- 
tron/ion mass ratio m e /mi [6|. When /? <C m e /rrii (i.e. 
when Alfven speed is much greater than electron and 
ion thermal speeds, Va S> v t h.i,vth.e) dominant mech- 



anism for sustaining E\\ is the parallel electron inertia 
and such waves are called Inertial Alfven Waves. When 
(3 S> m e /mi, (i.e. when Va <C vth,i,v t h,e) then thermal 
effects become important and the dominant mechanism 
for supporting En is the parallel electron pressure gradi- 
ent. Such waves are called Kinetic Alfven Waves. 

The context of this study is related to the theoretical 
plasma physics processes operating in the particle accel- 
eration in solar flares. Ref.pTj's Introduction gives a good 
overview of both theoretical and observational unresolved 
issues. If the solar flare particle acceleration happens in 
the corona, in order to explain the observed X-ray flux as 
the electrons smash into the dense layers of the Sun, im- 
plausibly large acceleration region volumes and/or den- 
sities are required. Different ideas have been put for- 
ward: (i) Substantial re-acceleration in the chromosphere 
of electrons accelerated in and injected from the corona 
can greatly reduce the density and number of fast elec- 
trons needed to produce a X-ray flux pj; However, this 
still seems problematic because recent analysis of the evo- 
lution of a radio spectrum from a dense flare § (Bastian 
et al. 2007) shows that a significant fraction of the en- 
ergy in the energetic electrons can be deposited into the 
coronal loop, as opposed to into the chromosphere. The 
estimates show that the energy deposited in the corona 
can approach about 30% flare energy, (ii) The electric 
circuit formed by precipitating and returning electrons 
9, 10]. The effect was found from the Ampere law which 
handles the circuit of injected and returning electrons 
when simulated with Particle- In-Cell and Fokkcr-Planck 
simulations, (iii) When wave-particle interactions in non- 
uniform plasma are taken into account, the evolution 
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of the Langmuir wave spectrum towards smaller wave- 
numbers, leads to an effective acceleration of electrons. 
Thus, the time-integrated spectrum of non-thermal elec- 
trons shows an increase in super-thermal electrons, be- 
cause of their acceleration by the Langmuir waves [Tlj | . 
(iv) Solar flare triggered DAWs, as opposed to the elec- 
tron beams, propagating towards the solar coronal loop 
foot-points, and accelerating electrons along the propa- 
gation paths [HI]. Key ingredient in such approach is 
the transverse density inhomogeneity (i.e. edges of solar 
coronal loops), which enable DAW to acquire non-zero 
perpendicular to the magnetic field wavelength, compa- 
rable to the above mentioned kinetic scales. This results 
in the generation of parallel electric fields, that can effec- 
tively accelerate electrons flsWIEj ]. 



Introduction section of Ref. [15| gives an overview of 
the previous work on this topic in some detail. Here we 
mention a latest addition, Ref. which studies the in- 
teraction of an isolated Alfven wave packet with a plasma 
density cavity. Ref.[15J considered particle acceleration 
by DAWs in the transversely inhomogeneous plasma via 
full kinetic simulation particularly focusing on the effect 
of polarisation of the waves and different regimes (iner- 
tial and kinetic). In particular, Ref.ff5j studied particle 
acceleration by the low frequency (w = 0.3u C i) DAWs, 
similar to considered in Ref. |13l If4fl , in 2. 5D geometry, fo- 
cusing on the effect of the wave polarisation, left- (ion cy- 
clotron branch) and right- (whistler branch) circular and 
elliptical, in the different regimes inertial (/3 < m e /mi) 
and kinetic (/? > m e /mi). A number of important con- 
clusions were drawn, including (i) The fraction of accel- 
erated electrons (along the magnetic field), in the density 
gradient regions is 20%-35% in 2.5D geometry, (ii) While 
keeping the power of injected DAWs the same in all con- 
sidered numerical simulation runs, in the case of right 
circular, left and right elliptical polarisation DAWs with 
E y /E z = 6 (with x being the direction of the uniform 
background magnetic field) produce more pronounced 
parallel electron beams, (iii) The parallel electric field 
for solar flaring plasma parameters exceeds Dreicer elec- 
tric field by eight orders of magnitude, (iv) Electron 
beam velocity has the phase velocity of the DAW. This 
can be understood by Landau damping of DAWs. The 
mechanism can readily provide electrons with few tens of 
keV. (v) When in 2.5D case the mass ratio was increased 
from rrii/m e = 16 to 73.44, the fraction of accelerated 
electrons has increased from 20% to 30-35% (depending 
on DAW polarisation). This is because the velocity of 
the beam has shifted to lower velocity. As there are al- 
ways more electrons with a smaller velocity than higher 
velocity in the Maxwellian distribution, one can conjec- 
ture that for the mass ratio mi/m e = 1836 the fraction 
of accelerated electrons would be even higher than 35%. 

In the present work we focus on the 3D effects on par- 
ticle acceleration and parallel electric field generation. 
In particular, instead of ID transverse, to the magnetic 
field, density (and temperature) inhomogeneity, we con- 
sider the 2D transverse density (and temperature) inho- 



mogeneity in a form of a circular cross-section cylinder, in 
which density (and temperature) varies smoothly across 
the uniform magnetic field that fills entire simulation do- 
main. Such structure mimics a solar coronal loop which 
is kept in total pressure balance. Section II describes 
the model for the numerical simulation, while the results 
are presented in section III. We close the paper with the 
conclusions in section IV. 



II. THE MODEL 

The general observational context of this work in out- 
lined in Fig. I, which shows that a solar flare at the 
solar coronal loop apex triggers high frequency, 0.3w C i, 
DAWs which then propagate towards loop footpoints. 
We only consider a single DAW harmonic, however in 
reality flare may trigger a wide spectrum of waves with 
different wave- numbers and frequencies, prescribed by a 
turbulent cascade. Some examples (albeit longitudinal 
(Langmuir), not transverse (Alfvenic) typ e) of creating of 
such turbulence can be found in Ref. [17| . We conjecture 
that if the electron beam has non-zero velocity compo- 
nent, transverse to the background magnetic field, kinetic 
plasma instabilities, such as electron cyclotron maser will 
also generate the transverse (Alfvenic) turbulent cascade. 
Another issue to consider in the context of DAW excita- 
tion is the possibility of excitation of KAWs or IAWs 
by means of magnetic field-aligned currents - essentially 
electron beams drifting with respect to stationary ions 
[HI or fast ion beam excitation [f9j]. Ref.[20| considered 
the situation when KAWs are excited by current (fluid) 
instability. The instability condition for this excitation 
by current is satisfied even for small drift velocity, e.g 
V£> = O.Iva, when KAWs can effectively grow. How- 
ever, Ref. [20] did not include the resonant excitation of 
DAWs by the inverse Landau damping because its insta- 
bility condition requires a larger drift velocity, in general, 
larger than the Alfven velocity. In a different regime, the 
importance of the Landau (Cerenkov) resonance for the 
particle acceleration and parallel electric field generation 
by the DAWs was stressed by Refs. flR l2l| . 

In our model (see Fig. 1) the transverse density (and 
temperature) inhomogeneity scale is of the order of 30 
Debye length (Ad) that for the considered mass ratio 
rrii/m e — 16 corresponds to 0.75 ion inertial length c/ui p i. 
Existence of such thin loop threads, tens of cm wide, in 
the solar corona is of course open to debate. The finest 
loops observed with TRACE have a width of the order 
of 1000-2000 km and have a monolithic structure with 
a single temperature at any cross-section, thus presently 
not supporting the thin thread concept. The TRACE 
CCD camera has 0.5 arcsec pixels which is 366 km on 
the sun. Therefore it is plausible that the smallest ob- 
served 1 000 km wide monolithic structures are probably 
too close to the resolution limit. The future high spatial 
resolution space missions such as Solar Probe Plus may 
possibly reveal the loop sub-structuring. 



We use EPOCH (Extendable Open PIC Collaboration) 
a multi-dimensional, fully electromagnetic, relativistic 
particle-in-cell code which was developed and is used 
by Engineering and Physical Sciences Research Council 
(EPSRC)-fundcd Collaborative Computational Plasma 
Physics (CCPP) consortium of 30 UK researchers. We 
use 3D version of the EPOCH code. The relativistic 
equations of motion are solved for each individual plasma 
particle. The code also solves Maxwell's equations, with 
self-consistent currents, using the full component set of 
EM fields E x ,E y ,E z and B x ,B y ,B z . EPOCH uses SI 
units. For the graphical presentation of the results, here 
we use the following normalisation: Distance and time 
are normalised to c/u pe and ui p e, while electric and mag- 
netic fields to aj pe cm e /e and Lo pe m e /e respectively. When 
visualising the normalised results we use no = 10 16 m -3 
in the least dense parts of the domain (y = 0, y = y m , 
z = 0, z = z m ), which are located at the edges of the 
simulation domain (i.e. fix uj pe = 5.64 x 10 9 Hz radian 
on the domain edges). Here uj pe = n e e 2 /(eom e ) is the 
electron plasma frequency, n a is the number density of 
species a and all other symbols have their usual meaning. 
The spatial dimension of the simulation box is fixed at 
x = 5000 and y = z = 200 grid points for the mass ratio 
mi/m e = 16. This mass ratio value corresponds to the 
in the inertial Alfven wave (IAW) regime because plasma 
beta in this study is fixed at j3 = 2(vth,i / c) 2 (u> P i / 'uj C i) 2 = 
n (0,0)k B T/(B 2 /(2n )) = 0.02. Thus (3 = 0.02 < 
m e /mi — 1/16 = 0.0625; This is the maximal value that 
can be considered with the available computational re- 
sources. Using 720 processor cores, the presented in this 
paper numerical run took 9 days and 14 hours. The grid 
unit size is Xp. Here A^ = Dt/i.e/^pe = 5.345 x 10 -3 m is 
the Debye length (v t h, e — \J k B T '/m e is electron thermal 
speed). This makes the spatial simulation domain size 
of x = [0,x m ] = [0,26.727m], y = [0,y m ] = [0,1.069m] 
z = [0,z m ] = [0,1.069m]. Particle velocity space is re- 
solved (i.e distribution functions produced in V x ,V y ,V z 
directions) with 30000 grid points with particle momenta 
in the range +6.593 x 10~ 22 kg m s _1 . We use 2.4 x 10 9 
electrons and the same number of protons in the simu- 
lation. In principle, plasma number density n(y, z) (and 
hence u pe (y,z)) can be regarded as arbitrary, because 
we use io pe (0,0) in our normalisations. We impose con- 
stant background magnetic field Bq x = 320.753 Gauss 
along x-axis. This corresponds to Bq x = 1.0(cu pe m e /e), 
so that with the normalisation used for the visualisa- 
tion purposes, normalised background magnetic field is 
unity. This sets uj ce /uj pe = 1.00. Electron and ion 
temperature at the simulation box edge is also fixed at 
T(0,0,0) = T e (0,0,0) = Ti(0,0,0) = 6 x 10 7 K. This in 
conjunction with no(0, 0) = 10 16 m~ 3 makes plasma pa- 
rameters similar to that of a dense flaring loops in the 
solar corona. 



We consider a transverse to the background magnetic 
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FIG. 1: A conceptual sketch of the model. A solar flare 
launches DAWs from the solar coronal loop apex which rush 
down towards the photospheric footpoints, accelerating elec- 
trons along the propagation direction and heating ions in the 
transverse direction. The transverse inhomogeneity scale is 
~ 30 Debye length (Ad) that for the considered mass ratio 
rrii/me = 16 corresponds to 0.75c/ 'u} pi . 




FIG. 2: Shaded surface plot of n e (x = x m /2, y, z) at t = 0. 



field variation of number density as following 

n t (y) = n e (y) 
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Eq.(l) implies that in the central region (across the y 
and z directions), the density is smoothly enhanced by 
a factor of 4, and there are the strongest density gradi- 
ents having a width of about 30 Ad around the points 
r s = y/ y 2 + z 2 = 51.5Ad and r s 
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magnetic field, thus the simulation is weakly non-linear. 
Such transverse electric field E± — E y%z driving gener- 
ates L-circularly polarised DAW. For further details on 
polarisation see paragraph before Eq.(7) from Ref. |15j . 
In the case of L-circularly polarised DAW the dispersion 
is 



1/2 



{u ce + u){u> ci - bj) 



(5) 



Eq.(5) shows that when to = uj c 
nance occurs. 



the ion-cyclotron reso- 



FIG. 3: Shaded surface plot of T e (x = x m /2, y, z) oc l/n e (x — 
x m /2, y, z) at t = 0. 

r = (j/o, zo) = (IOOAd, IOOAd). Fig. 2 shows n e (x = 
x m /2, y, z) at t = 0, which indicates that the number 
density at the cylindrical tube centre increases by factor 
of 4 compared to the exterior. This density behaviour 
represents the solar coronal loop. The background tem- 
perature of ions and electrons are varied accordingly 

Ti(y, z)/T = T e (y, z)/T = f(y, z)' 1 , (2) 

such that the thermal pressure remains constant. Be- 
cause the background magnetic field along the x- 
coordinate is constant, the total pressure is also constant. 
Fig. 3 shows T e (x = x m /2,y,z) oc l/n e (x = x m /2,y,z) 
at t = 0. Note how the temperature at the cylindrical 
tube centre decreases by factor of 4 compared to the ex- 
terior. 

The DAW is launched by driving domain left edge, 
x = as follows 

E y (l,y,z,t + At) =E y (l,y,z,t)- 
A y sm{uj d t) (1 - cxp [-(i/i ) 2 ]) , (3) 

E z (l,y,z,t+ At) = E z (l,y,z,t)- 
A z cos{u d t) (1 - exp [-(<Ao) 2 ]) , (4) 

which produces the case of left-polarised DAW, i.e. ion 
cyclotron wave. u>d denotes the driving frequency which 
was fixed at u>d = 0.3w C i. The fact that ujd < <jJ C i ensures 
that no ion-cyclotron damping is present and also that 
the generated DAW is well on the Alfven wave branch 
with dispersion properties similar to the Alfven wave (cf. 
Fig. 3.4 from Ref. [22j] ) . to is the onset time of the driver, 
which was fixed at 3.0W" 1 i.e. 48.000w~ e 1 for the case 
of nii/m e — 16. Thus the driver onset time is (9 
comes from 3 2 ). The initial amplitudes of the E± are 
chosen as A y = A z = 0.05(cBq x ). Therefore magnetic 
perturbation amplitude of DAW is 5% of the background 



III. RESULTS 

Fig. 4 presents the numerical simulation results for the 
electromagnetic fields. We see from Figs. 4(e) and 4(f) 
that in response to the transverse electric field driving, 
which proceeds according to Eqs.(3) and (4), left-hand, 
circularly polarised IAW has developed three wavelength. 
If we use Eq.(5), then for the considered plasma pa- 
rameters, 3 x (2ir/k)/(c/u! pe ) = 201.7. We gather from 
Fig. 4(e) that the waveform at y = has travelled to 
x pa 200. It has not reached z — 201.7 quite yet be- 
cause the driving onset time is 9W" 1 . We see that E y 
component develops phase mixing i.e. waveform dis- 
torts across the transverse y and z coordinates generat- 
ing strong transverse gradients. The transverse magnetic 
fields have similar behaviour to transverse electric fields. 
However, phase- mixed field component is B z (Fig. 4(c)). 
Movie that shows time dynamics of B z (x,y,z — z m /2) 
(a plane that cuts through the middle of the over-dense 
tube, along the background magnetic field) is shown in 
the supplementary material [23[. Fig. 4(d) shows that 
parallel electric field is generated only in the density in- 
homogeneity regions 4 < r < 7 and 13 < r < 16, where 
r = \Jy 2 + z 2 . The time evolution of this parallel electric 
field is shown in two different cross-sections in the supple- 
mentary material [23| : E x (x,y,z = z m /2) (a plane that 
cuts through the middle of the over-dense tube, along the 
background magnetic field) and E x {x — x m /8,y, z) (a 
plane that cuts through the 1 / 8th of the over-dense tube, 
across the background magnetic field). We gather from 
movie 2 [23[ that in this 3D simulation, in the plane cut 
along the background magnetic field, behaviour is similar 
to 2.5D results both in PIC simulation (see e.g. Fig. 1 
from Ref. [l5j ) or cold plasma approximation (see anima- 
tion 2 (in Supplementary Data) from Ref. 24] ) . What is 
new in 3D is the ability to study the transverse structure 
of the parallel electric field. Somewhat unexpected re- 
sults are shown in movie 3 [23j . The bright yellow blob 
that can be seen rotating around the centre is a marker 
to guide ones eye. It rotates in the clockwise direction 
with the angular frequency of 0.3aj c ;. Its dynamics is 
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FIG. 4: Contour (intensity) plots of the following physical quantities at t en d = 75w ci : (a) B x (x,y,z = z m /2) — Bo, (b) 
B y (x,y,z = z m /2), (c) B z (x,y,z = z m /2), (d) E x (x,y,z = z m /2), (e) E y (x,y,z = z m /2), (f) E z (x,y,z = z m /2). These 
correspond to a vertical cut through the mid-plane 2 = z m /2. 



described by the equation 

marker (y, z,t) = 2.0 x 10~ 2 exp^ 

-{[(!/ - J/m/2) COs(0.3w CJ i + 7r)- 

(z - Zm/2) sin(0.3w ci i + tt) - y m /4]/0.5} 2 
— {[(z - z m /2) cos(0.3w c rf + 7r) + 

(2/ - y m /2) sin(0.3w C2 i + tt) - z m /4]/0.5} 2 ) . (6) 

We see only the marker blob rotating in the clockwise 
direction until about time t < 600w pe . Then the gener- 
ated parallel electric field reaches the plane correspond- 
ing to x — x m /8 in which this movie is produced. From 
t > 600cj pe until the end of the simulation we see that 
parallel electric field is localised within circular bands 
4 < r < 7 and 13 < r < 16 and rotating in the clock- 
wise direction with the frequency 0.3w C i. This seems 
initially surprising why the parallel electric field should 
rotate. However, if one recalls that the IAW also ro- 
tates in the same clockwise direction, then this seems 
logical (that both parts of the wave E± that is used as 
a driver and the generated E\\ both rotate in the same 
sense and with same angular frequency) . Note that rota- 
tion in movie 3 and translatory motion in movies 1 and 
2 [23| appear step-like because there are only 20 time 
snapshots taken throughout the simulation that ends 



at tend — T^w^ 1 = 1200^^. Thus in this animation 
At = 3.750;^ . Also, in all movies distances on x — y 
axes are quoted in number of grid points - not in c/u> pe 
as in the rest of the paper. Thus, the conclusion here 
is that the generated parallel electric field is localised in 
the density inhomogeneity region and rotates in the same 
direction and with the same angular frequency as the ini- 
tially launched DAW. 

We quantify the particle acceleration and plasma heat- 
ing by defining the following indexes: 

AH te (t) = 

J'\Z\>(v th ,.) fe(vx,t)dv x / (x m ir[(r - r ) 2 - (n - r ) 2 ]) 

I-00 fe( v x, 0)dv x / (x m y m Z m ) 

AH XA (t) = 

■f\7 ± \>(v th:i ) /»(«!-> *) d «x/ {x m -K[{r - r ) 2 - (n - r ) 2 ]) 

I-00 f* ( w± ' / (x m y m z m ) 

where r Q and are outer and inner radii of the density in- 
homogeneity shown in Fig. 2 and rg is its mid-point (cen- 
tre). f e< i are electron or ion velocity distribution func- 
tions. They are normalised in such a way that when in- 
tegrated by all spatial and velocity components they give 
the total number of real electrons and ions (not pseudo- 
particles that usually mimic much larger number of real 
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FIG. 5: Time evolution (at t = and t — t en d = 75W" 1 ) of electron and ion velocity distribution functions versus velocity x, y and 
z components on a log-linear plot: (a) f e (v x ,t = 0) dashed (inner) curve and f e (v x , t = t en d) solid (outer) curve, (b) f e (v y ,t = 0) 
dashed (inner) curve and fe(v y ,t — tend) solid (outer) curve, (c) f e (v z ,t = 0) dashed (inner) curve and f e (v z ,t = t en d) solid 
(outer) curve, (d) fi(v x ,t = 0) dashed (inner) curve and fi(v x ,t = t en d) solid (outer) curve, (e) fi{v y ,t = 0) dashed (inner) 
curve and fi{v y ,t = t en d) solid (outer) curve, (f) fi(v z , t = 0) dashed (inner) curve and fi(v z ,t — t ent j) solid (outer) curve. Time 
evolution (at 20 time intervals between t — and t = t en d) of the following: (g) max \E x (x, y, z, t)|, triangles connected with 
a solid curve and (h) AH\^ e (t) index, diamonds connected with dashed curve, according to Eq.(7), AHx,i(t) index, triangles 
connected with a solid curve, according to Eq.(8). 



particles in the PIC simulation). The definition used in 
Eq.(7) effectively provides the fraction (the percentage) 
of accelerated electrons, parallel to the magnetic field, 
in the density gradient regions. Eq.(8) gives the frac- 
tion (the percentage) of heated ions, perpendicular to the 
magnetic field. In the case of 2.5D similar indexes have 
been defined in Ref. 15]. We refer interested reader for 
an underlying discussion of these indexes in the former 
reference. We gather from Fig. 5(a)-5(c) that electron 
acceleration is confined only to the parallel to the back- 
ground magnetic field direction. This can be evidenced 
by appearance of the bumps in the f e (v x ) at t = t en d, 
which correspond to the electron beams produced by the 
generated E x — E\\ in the density inhomogeneity. By 
comparing Fig. 5(a) (electrons) with Fig. 5(d) (ions) we 
see that no ion beams are generated. Instead we see from 
Figs. 5(e) and 5(f) that ions are heated in the perpen- 
dicular direction as the distribution function at t = t enc [ 
widens (whilst no ion beams are produced). Note that 
such transverse ion heating of plasma has been studied 
in the past too [25j-|28|, while the main focus here is the 
electron acceleration in the parallel direction. Compar- 



ing Fig. 5(a) (electrons) with Fig. 5(d) (ions) also gives us 
an indication that because in the density gradient region 
(circular ring) ions essentially do not feel IAW driving 
in the parallel direction while electrons are accelerated, 
effective charge separation takes place. This is consis- 
tent with the cold plasma approximation, see animation 3 
(in Supplementary Data) from Ref.[13|. We gather from 
Fig. 5(g) that the value of the generated parallel electric 
field, in normalised units, is 0.07. It can be also seen 
from Fig. 5(h) that AH» e (t) and AH± y i(t) indexes attain 
values of 0.45 (45%) and 0.75 (75%) respectively. Fig. 3 
from Ref.yJI directly corresponds to our Fig. 5 here but 
now for the 3D case. Thus we can make some compar- 
isons. The differences are: (i) in 3D case parallel electron 
beams are more pronounced, as the solid line in Fig. 5(a) 
attains higher values also over the larger v x range, (ii) in 
3D case the transverse ion heating is more pronounced 
too, as the solid lines in Figs. 5(e) and 5(f) reach higher 
values compared to the 2.5D case, (hi) in 3D case parallel 
electric field attains values roughly a factor of two larger, 
as Fig. 5(g) shows that the value of the parallel electric 
field attains 0.07 versus 0.03 in the 2.5D case, (iv) the 
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FIG. 6: Time evolution (at t = and t — t en d = 75uj~ 1 ) of electron and ion velocity distribution functions (with positive 
velocities) versus kinetic energy of particles £ on a log- log plot. Top row: (a) f e ,Vm(E,t = 0) thin solid (inner) curve and 
f e ,Vx(E,t = t en d) thick solid (outer) curve, (b) fe,Vy{E,t = 0) thin solid (inner) curve and f e ,Vy(E,t — tend) thick solid 
(outer) curve, (c) j e ,Vz{E,t = 0) thin solid (inner) curve and f e ,Vz(E,t = tend) thick solid (outer) curve. Vertical line in 
(a) corresponds to the kinetic energy, Edaw = (7 — l)m e c 2 , with 7 = 1/Vl — 0.201 2 and 0.201c being phase velocity of 
left-hand polarised IAW in the low density region. Dotted line is 2 x W 10 E~ 15 and dashed line is 10 7 -E~ 4 . Bottom row: (d) 
fi,Vx(E,t — 0) thin solid (inner) curve and fi t vx{E,t = tend) thick solid (outer) curve, (e) fi,v y (E,t — 0) thin solid (inner) 
curve and fi,v y (E,t = t en d) thick solid (outer) curve, (f) fi,vz(E,t — 0) thin solid (inner) curve and fi,vz(E,t = t en d) thick 
solid (outer) curve. 



fraction of accelerated electrons, parallel to the magnetic 
field, in the density inhomogeneity is factor of two larger, 
as it can be seen from Fig. 5(h) that AH» e (t en d) — 0.45 
versus AHn e (i e nd) = 0-2 in the 2.5D case and the frac- 
tion of heated ions, perpendicular to the magnetic field, is 
also factor of two larger, as it can be seen from Fig. 5(h) 
that AH± ii (t en( i) — 0.75 versus AHj_^(t end ) = 0.35 in 
the 2.5D case. 

Fig. 6 presents the same data in panels (a)-(f) from 
Fig. 5 but now in terms of kinetic energy E = (7 — l)m e c 2 
and also plotted on the log-log plot. Note that in Fig. 5 
log-normal plot was used. To produce Fig. 6 only v x > 
branch was used. The reason for producing Fig. 6 is two- 
fold: (i) to demonstrate that the electron spectrum in 
the parallel to the field direction shows so-called "knee" 
- a break in the spectra power-laws; (ii) to demon- 
strate that the "knee" can be interpreted is the Landau 
damping (Cerenkov resonance effect) due to the wave- 
particle interactions. There are several explanations for 
the observed double power-laws in the electron spectrum. 
By double we mean one (usually smaller oc E~ a with 
2 < a < 4) power-law index below a certain fixed en- 



ergy E < E br (where 50 keV < E br < 150 keV and then 
another (usually larger oc E~ a with 4 < a < 7) power- 
law index for E > E br 29]. Chapter 13.3.4 of Ref.[30| 
quotes several explanations for the double power-law also 
known as the "knee". The most common explanation is 
collisional relaxation of an electron beam, with added 
refinements such as the generation of Langmuir waves 
and plasma density inhomogeneity [ljj . Here we offer 
an alternative explanation in that flare generated DAWs 
produce electron beams on the transverse density inho- 
mogeneities (magnetic loop edges) which in turn produce 
spectra with the structure that resembles the knee. 

Fig. 7 is analogous to Fig. 6 but it has been produced 
with the simulation data from run L16 from Ref.[l5j that 
corresponds to 2.5D. Comparison of Fig. 7 (2.5D case, 
Ref . [la] ) and Fig. 6 (3D case, this paper) enables to grasp 
the differences between 2.5D and 3D geometry. It is evi- 
dent from Figs. 7 and 6 that in the 3D case acceleration 
of electrons along the magnetic field and heating of ions 
in the transverse direction is more efficient. Also, the 
knee feature in the parallel to the magnetic field electron 
spectra can be seen. 
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FIG. 7: As in Fig. 6 but for run L16 from Ref.QJ]. Dotted line is 10 9 £T 2 and dashed line is 2 x 10 3 £T 6 . 




FIG. 8: As in Fig. 6 but for run L73 from Ref . [l5l] . Vertical line in (a) corresponds to the kinetic energy, Eoaw — (7~ l)m e c 2 , 
with 7 = — 0.097 2 and 0.097c being phase velocity of left-hand polarised IAW in the low density region. Dotted line is 

10 10 iT 2 and dashed line is lO^B -8 . 
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Fig. 8 is analogous to Fig. 6 but it has been produced 
with the simulation data from run L73 from Ref.[l5j that 
corresponds to 2.5D. Comparison of Fig. 8 (2.5D case, 



Ref. 15 1 mass ratio mi/m e — 73.44), Fig. 7 (2.5D case, 
Ref. 1151] mass ratio mi/m e — 16) and Fig. 6 (3D case, 



this paper, mass ratio rrii/m e = 16) enables to gauge the 
differences between 2.5D versus 3D geometry as well as 
one mass ratio versus another (note that considered mass 
ratios land the physical system into different regimes. 
Case of mi/m e — 16 corresponds to IAW regime, whereas 
mi/m e — 73.44 corresponds to KAW regime). We gather 
from Fig. 8 that double power-law fit to the parallel elec- 
tron spectrum is better compared to Figs. 7 and 6, possi- 
bly indicating that the KAW regime provides more dou- 
ble power-law-like behaviour. It should be noted that in 
Fig. 7 and 6 corresponding to IAW regime power-law fit 
seems somewhat arbitrary. This partly could also be due 
to the fact that the observations usually provide time- 
averaged electron spectrum whereas no-time averaging 
has been used in Figs. 6-8. A more profound result is 
that the energy, Ed aw, that corresponds to the break 
in the power-law, E^, approximately corresponds the en- 
ergy of an electron which moves with the phase speed of 
the DAW, V p h- In other words we note that as we in- 
crease the mass ratio from 16 to 73.44 both the break 
in the spectrum and V p h decrease. This demonstrates 
that the "knee" can possibly be interpreted as the Lan- 
dau damping (Cerenkov resonance effect) of DAWs due 
to the wave-particle interactions. 



IV. CONCLUSIONS 

The aim of this work was to explore the novelties 
brought about by 3D geometry effects into the problem of 
particle acceleration by DAWs in the solar flare and also 
to Earth magnetosphere auroral zone studied in an earlier 
2.5D geometry work [151 ] . It should be noted that here we 
have studied the case of plasma over-density, transverse 
to the background magnetic field, whereas references that 
deal with Earth auroral zone usually consider the case of 
plasma under-density (a cavity), as dictated by the dif- 
ferent applications considered (solar coronal loops and 
Earth auroral plasma cavities). As far as the genera- 
tion of parallel electric field and associated particle ac- 
celeration is concerned, there is no difference whether the 
transverse density gradient is positive or negative - for 
the mechanism to work it has to be non-zero. This is 



because whistlers (and ion cyclotron waves) that prop- 
agate strictly along the magnetic field display no Lan- 
dau damping, since the longitudinal component of the 
wave electric field is zero [31(, page 124. The latter be- 
comes non-zero as the wave front turns due to the trans- 
verse density inhomogeneity. Thus we investigated a pro- 
cess of particle acceleration by left-hand, circularly po- 
larised inertial Alfven wave in a transversely inhomoge- 
neous plasma, using 3D particle-in-cell simulation. We 
considered a cylindrical tube that contains a transverse 
to the background magnetic field inhomogeneity with a 
scale of the order of ion inertial length. Such numerical 
setup is relevant for solar flaring loops and earth auroral 
zone. In such structure IAWs with frequency 0.3w C i were 
launched and allowed to develop three wavelength. The 
following key points have been established: Propagation 
of IAW in such a system generates time-varying paral- 
lel electric field, localised in the density gradient regions, 
which accelerate electrons in the parallel to magnetic field 
direction. Perpendicular electric field of IAW also effec- 
tively heats ions in the transverse direction. The gener- 
ated parallel electric field rotates in the same direction 
and frequency as the " parent" IAW. The fully 3D kinetic 
simulation demonstrates electron acceleration efficiency 
in the density inhomogeneity regions, along the magnetic 
field, is of the order of 45% and ion heating, in the trans- 
verse to the magnetic field direction, is about 75%. The 
latter is a factor of two times higher than the previous 
2.5D analogous study and is broadly in agreement with 
the solar flare particle acceleration observations. Log- 
log plots of electron spectra seem to indicate that the 
"knee", frequently seen in the solar flare observations, 
can be interpreted is the Landau damping of IAWs due 
to the wave-particle interactions. 
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